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Determination of the position of a pulse peak 
FIELD OF THE INVENTION 

The invention relates to a method for determining the 
position of the peak of a pulse in a signal received at a 
receiver. The invention relates equally to a device and 
to a cellular communication system which realize this 
method. 

BACKGROUND OF THE INVENTION 

The position of the peak of a pulse in a received signal 
can be required for example for determining the delay of 
the signal when propagating from a transmitting unit to 
the receiver. 

The delay of signals at a receiver can be evaluated for 
instance by a location service for determining the 
current location of the receiver. In case of a line-of- 
sight transmission, the delay of a signal is directly 
dependent on the distance between the receiver and the 
respective transmitting unit. 

Such a location service can be provided in particular by 
a CDMA (code division multiple access) based satellite 
positioning system or by a CDMA, based terrestrial 
cellular positioning system. In a CDMA based system, a 
data sequence is used by a transmitting unit to modulate 
a sinusoidal carrier, and then the bandwidth of the 
resulting signal is spread to a larger value, e.g. by 
multiplying the modulated signal with pseudo-random bits 
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*riv*d from a ^ epreading codSi 
usually r.f erred to as chips. 

dentin! 71T ^ SeWhin9 PrOCStee *~ f °^ «* 
detecting a delayed slgnal taking ^ ^ 

Propagation path is normally carried out in the imn„i 

units The delay can he estimated e.g. by an edge 
detects m the impulse response J 
"gnale. The length of the impulse response profile i. ^ 
much longer than the „i dt h of the signal shape 
TZT° Z " ^ " arCh ^ Sd9e 1S P^°-ed along the 
signal b r ePOnae : ""^ ^ 3 ""^ ~*«- °* «- 
pra-defaned threshold. The edge detection is thus a 
hitting process. The threshold has to he est on the one 
hand high enough in order to avoid that a noise peak is 
detected a S signal edge, which would result in a false 

e'nouTh T ^ " her hand ' thr " h0ld 1135 " b * "* 1- 

anough to guarantee that the signal edge is detected aven 

if the signal strength is rather weak. 

«• delay of a signal can only be determined accurately, 
"hen the exact ^sition of the peak of a signal pulse L 
known since this is the only claar reference point in 
the pulse. The result of the hitting process, however, is 
usually a pcsition on the left side of the signal peak, 
i.e. on the side which is closer to a delay of 
since a pulse will usually be detected before its peak is 
reached. Thus, the error of the edge detection is 
negatively biased. The error is mors related to the 
signal than to the SN * (signal-to-noise ratio, , which 
mean, that the variance can be high. The hit will occur 
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between close to the peak of the pulse for weak signals 
and close to the bottom of the pulse for strong signals. 

In case of a sampling rate of 2 samples per chip, the 
delay estimation error can therefore range from 0.0 chips 
to -i.o chips. The average error is then about -0.5 chip 
for a triangular shape or waveform of the pulse. 

» 

Since the error has a negative bias, the simplest way to 
reduce the error is to introduce a positive factor to 
compensate for the bias. If the signal level is higher, 
the error is also bigger. Therefore, the compensation 
factor should be adaptive to the signal level. Still, 
with such a general compensation, a significant average 
error remains. 

SUMMARY OF THE INVENTION 

It is an object of the invention to improve the accuracy 
of the determined position of the peak of a pulse. 

This object is reached according to the invention with a 
method for determining the position of a peak of a pulse 
in a signal received at a receiver, which method 
comprises as a first step taking samples of the received 
signal. In a next step, at least three samples are 
determined, of which at least one has a signal strength 
exceeding a threshold value. The samples can, but do not 
have to be consecutive. Then, the position of the pulse 
peak is determined based on an interpolation of at least 
two of the determined samples. These at least two samples 
are selected based on the signal strengths of the at 
least three determined samples. It is further proposed 
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that the interpolation includes an evaluation of the 
signal strength of the at least two samples. 

The object of the invention is moreover reached with a 
device comprising means for realizing the proposed 
method. The device can be for example the receiver, i.e. 
the receiver itself computes the position of received 
pulse peaks. Alternatively, the device can be external to 
the receiver. In this case, the receiver has to provide 
information on received signals to this device, in order 
to enable the device to perform the proposed processing. 
In case the receiver is a mobile terminal or integrated 
in a mobile terminal, such an external device can be in 
particular a network element of a cellular communication 
network. Finally, the object of the invention is reached 
with a cellular communication system comprising the 
proposed device. 

The invention proceeds from the idea that the position of 
a pulse peak can be determined more accurately, if the 
position is not determined simply by an edge detection 
searching for a single sample which exceeds a threshold 
value, but by an interpolation of at least two samples, 
of which at least one exceeds a threshold value. The 
proposed interpolation evaluates the signal strengths of 
the samples used for the interpolation, thereby taking 
account of the difference in the distance of the samples 
from the peak of the pulse. In addition, the at least two 
samples can be selected in most cases such that they can 
be assumed to lie advantageously on opposite sides of the 
pulse peak, if the signal strength from an appropriate 
number of samples is taken into account , 
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It is an advantage of the invention that the estimation 
accuracy of the position of the peak is increased. 
Thereby, e.g. the delay of received signals can be 
estimated more accurately. 

Preferred embodiments of the invention become apparent 
from the dependent claims . 

Advantageously, the equations for the interpolation are 
selected based on the employed model of the pulse shape. 
For example, in case the model of the" pulse shape is a 
triangle, a linear relation between the distance of the 
samples to the peak position and the signal strengths of 
the samples can be assumed for the interpolation. 

A model of the pulse shape which has a triangular shape 
results in a particularly simple approach. Even though in 
practice the real pulse shape will usually not correspond 
to a triangle, the triangle approach leads to quite 
accurate peak position estimates for other band-limited 
signal pulse shapes as well. While different pulse shapes 
have a different behavior in the delay estimation, the 
interpolation itself is a shape-fitting procedure. 

In order to compensate in addition for* the differences 
between the selected model and the real pulse shape, 
either the samples used in the interpolation can be 
weighted with different weighting coefficients or the 
resulting estimate can be adjusted. The compensation can 
be based on known deviations between the model of the 
pulse shape and the real pulse shape and/or on the signal 
strength of samples. 
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Preferably, the method according to the invention 
differentiates between different situations that may- 
occur. Thereby, the optimal samples can be selected for 
the interpolation in each situation. With such a 
differentiation, the method according to the invention is 
also suited to deal with multipath propagation of 
signals. The differentiation can be based in particular 
on the differences in the successions of the signal 
strengths of the considered samples . 

The invention can be employed e.g. for supporting the 
acquisition of signals, for instance in any kind of 
location service. 

It can be used in particular, though not exclusively, for 
determining the delay of signals in a satellite 
positioning system or in a cellular positioning system. 
The satellite positioning system can be for instance GPS 

(global positioning system) or Galileo, while the 
cellular positioning system can be implemented for 
example in a 3G (3rd generation) communication system, 
like a GSM (global system for mobile communications) , 
CDMA and/or GPRS (general packet radio system) system. 

BRIEF DESCRIPTION OF THE FIGURES 

Other objects, features and advantages of the present 
invention will become apparent from the following 
detailed description considered in conjunction with the 
accompanying drawings, wherein 

Fig. 1 is a flow chart illustrating an embodiment of the 
method according to the invention/ 
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Fig. 2 shows an exemplary model of a pulse shape 
employed in the method of figure 1; 

Pig. 3 illustrates a first situation dealt with in the 
method of figure 1; 

Fig. 4 illustrates a second situation in the method of 
figure 1; 

Fig. 5 illustrates a third situation in the method of 
figure 1; and 

Fig. S shows an exemplary real pulse shape together with 
a model of a pulse shape. 

DETAILED DESCRIPTION OF THE INVENTION 

Figures 1 is a flow chart which illustrate a first 
embodiment of the method according to the invention, 
which is implemented in an exemplary CDMA system. The 
CDMA, system comprises a plurality of base stations of a 
cellular communication network and a mobile station which 
is able to communicate via the air interface with these 
base stations. The method is used for estimating at the 
mobile station the delay of signals received from some of 
the base stations, in order to determine the current 
position of the mobile station. 

For the method of figure 1, the pulse shape in the 
impulse response of a received CDMA signal, which has 
passed a matched filter, is supposed to be triangular and 
is thus also modeled to be triangular. The triangular 
model is depicted in figure 2. The triangle 21 is 
isosceles and has a duration of 2 chips, since the 
coverage of a peak in a CDMA system is 2 chips. 

In a first step in the method of figure 1, samples are 
taken at the mobile station from the impulse response 
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with a sampling rate of two samples per chip. Obviously, 
a sampling rate of more than two samples per chip could 
be used as well. In this case, however, the equations 
presented in the following would have to be adjusted 
accordingly. 

The taken samples are then compared to a predetermined 
threshold value, starting at a position which corresponds 
to a delay of zero. The threshold value is selected such 
that the false alarm rate is below a desired value, while 
it is ensured at the same time that a real peak is 
detected with a desired probability. The threshold value 
may be varied to this end based e.g. on the signal 
strength of received signals. 

r 

The first sample at a position x 0 which exceeds the 
threshold value is considered as a "hit n , i.e. as a 
detection of the edge of the signal pulse in the impulse 
response which took the shortest propagation path. 

In the following, a differentiation between three 
possible situations is carried out. 

In case the strength A(x 0 ) of the first sample at position 
x 0 is smaller than the. strength A(x x ) of the second sample 
at position x x , there are two possible situations , 

In the first possible situation Si, the strength A(xi) of 
the second signal at position Xi is larger than the 
strength A(X2) of the third signal at position x 2 . The 
first situation Si is thus give by the following set: 

5, €[4x 0 )<4x I )]nWx 1 )^4^)] 
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In the first situation, a single path has been detected, 
and the peak is located somewhere between position x 0 and 
position x 2 . In this case, the three values A(Xo) , A(x x ) 
and A(x 2 ) can be used for an interpolation to find a more 
accurate peak position. 

An example for this situation is illustrated in figure 3 . 
The figure is a diagram depicting the strength A of the 
three samples at the respective position x<>, Xi, x 2 of the 
samples- In addition, the triangular model 31 of the 
pulse shape is indicated with dotted lines. The signal 
strength of the three samples is equal to the value of 
the triangle at the corresponding position. It can be 
seen that the first sample is positioned on the left hand 
side of the peak 32, while the second and the third 
sample are located on the right hand side of the peak 32. 
In another example for this situation, the second sample 
could also be located on the left hand side of the peak. 
A zero delay is somewhere to the left of the triangle 31. 

Por the interpolation, the "early minus late" technique 
can be applied. With this technique, the estimated peak 
position Xp, and thus the delay, is given by: 



where the unit of second term is chips. 

In figure 3, the estimation of the peak position Xp 
resulting for the depicted example is indicated. If the 
signal is strong, a perfect interpolation for the delay 
estimation can be obtained. The error is SNR dependent. 
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In the second possible situation S 2 , in contrast , the 
strength A(xi) of the second signal at position x x is 
smaller than the strength A(x 2 ) of the third signal at 
position x 2 . The second situation 3 2 is thus give by the 
following set : 

S 2 [A(x Q )< A( Xl )]f] [A{x x ) < A(x 2 )] 

Such a distribution of sample strengths is indicative of 
a multipath propagation resulting in a so-called "fat 
finger" . The second situation results more specifically, 
when a pulse transmitted by a base station propagates via 
several paths to the mobile station, and the first pulse 
arriving at the mobile station, e.g. a pulse which 
propagated on a line -of -sight path, is followed in an 
overlapping manner by a second pulse. 

An example for the second situation is illustrated in 
figure 4, Figure 4 is again a diagram depicting the 
strength of the three samples at the respective position 
Xo, x a.' x 2 of the samples. Here, two overlapping triangles 
41, 43 f each modeling a pulse arriving at the mobile 
station, are indicated in addition with dotted lines. The 
first sample, having a signal strength of A(x 0 >, is 
positioned on the left hand side of the peak 42 of the 
first pulse 41 at position x 0 , and the second sample, 
having a signal strength of A(x x ) , is positioned on the 
left hand side of the peak 42 of the first pulse 41 at 
position x x . The strength A(x a ) of the subsequent third 
sample at position x 3 is larger than the strength A(xi) of 
the second sample, since it comprises a component of the 
second pulse 43. 



02/12 02 LUN 16:40 [N° TX/RX 6019] 



2. D EZ. 2002 16:43 COHAUSZ & FLORACK "NR. 850 "G. 10 7 

PCT/IB02/05059 



- 11 - 



Therefore, the value of the third sample at position x 2 is 
not reliable. Nevertheless, the first and the second 
sample can be used to reduce the inaccuracy of the first 
sample. The position x p of the peak 42 of the first pulse 
41 can be estimated by: 

where the unit of second term is chips. 

For this estimation, there are two error bonds given by 
the two extreme cases. In the first extreme case, the 
position xc of the first sample is near the top 42 of the 
first triangle 41, and the maximum error caused by the 
lack of information on a reliable third sample is +0.25 
chips. In the second extreme case, the position x 0 of the 
first sample is near the foot of the first triangle 41, 
and the maximum error caused by the lack of information 
on a reliable third sample is -0.25 chips. 

The multipath interference may cause an additional error 
of +0.25 chips in the maximum, when the hit is near the 
top of the first path signal and the following path 
signal is very strong, i.e. A(x x ) » A(x 0 ) . This multipath 
error in the delay estimation is always positive. 

In case the strength A(x 0 ) of the first signal at position 
xo is not smaller than the strength A(x^) of the second 
signal at position Xi, the third possible situation S 3 is 
given, which can be expressed by the following: 



*°4 
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5 3 eM*o)M*i)] 

In this third situation, the first sample at position x 0 
can be assumed to be close to the peak. Still, the first 
sample might not be positioned exactly at the position of 
the peak. 

An example for the third situation is illustrated in 
figure 5. The figure is a diagram depicting the strength 
of the three samples at their positions x 0 , x x , x 2 . In 
addition, a triangular model 51 of the pulse shape is 
indicated again with dotted lines. It can be seen that 
all three samples are positioned on the right hand side 
of the peak 52 of the triangle 51. The third sample at 
position x 2 did not result in a signal strength exceeding 
the predetermined threshold value. 

In order to obtain an information on the exact peak 
position Xp, it has to be ensured that information from 
both sides of the peak 52 is available. To this end, in 
addition an earlier sample at position x.i is taken into 
account. The signal strength A(x. a ) of the preceding 
sample at position x_i is indicated in figure 5 as well. 
The interpolation is then carried out between the earlier 
sample position x- x and the second sample at position Xi. 
Thereby, the situation corresponds basically to the first 
described situation. 

As in the first situation, the "early minus late" 
technique can be applied for the interpolation. The 
estimated peak position x p is then given by; 
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where the unit of second term is chips. 

In the example of figure 5, the resulting estimation for 
the peak position Xp lies between position x-i and 
position Xo- 

Summarized, the estimation error can be limited to a 
value between -0.25 and +0.25 chips with the proposed 
method in any situation. The multipath error may add with 
a maximum of +0.25 chips. There is no signal -strength 
dependent bias required for the estimation. The stronger 
the signal, the smaller the error. The averaged error can 
be much smaller than 0.25 chips for strong signals. This 
corresponds in the case of a cellular 3G network to an 
estimation accuracy of ±20m, if there is no multipath 
propagation, and of -2 0m to +40m in the presence of 
multipath propagation. 

In a second embodiment of the invention, again a 
triangular model of the pulse shape is used, but the 
actual pulse shape is known to be a square -root -raised 
cosine pulse shape, which can be expressed by the 
following equation: 



Figure 6 shows a corresponding square -root -raised cosine 
pulse shape 65 with a solid line, and the triangular 
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model 61 with a dashed line. It can be seen that the 
differences between the two shapes 61, 65 within ±1 chips 
are quite small. The error caused by the shape 
differences is less than 5% in this case. 

In order to further improve the delay estimation, 
deviations of the model of the pulse shape 61 from the 
real pulse shape 65 are taken into account in addition to 
the steps described with reference to figure 1 , 

More specifically/ the samples of which the signal 
strengths A(x) are is to be used in the interpolation are 
first weighted with different weighting coefficients 
compensating for the differences in shape , as will be 
explained in the following by way of example • 

The situation presented in figure 6 corresponds basically 
to the first situation presented in figure 3. Thus, a 
first sample at position x 0 and a third sample at position 
x 2 are used for the actual interpolation. The first sample 
at position x 0 has a measured signal strength A(x 0 ) which 
exceeds the value of the triangle 61 at this position. 
Thus, this strength A(Xo) is weighted lower. The third 
sample at position x 2 , in contrast, has a measured signal 
strength A(x 2 ) which is lower than the value of the 
triangle 61 at this position. Thus, this strength A(x 2 ) is 
weighted higher* As a result, the square-root raised 
cosine curve is approached to the triangle curve, and the 
interpolation principle for the triangle can be also 
applied to the square-root raised cosine pulse shape. 

The interpolation is then carried out as described above 
with reference to figure 1 with the weighted samples. 
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In a third embodiment of the method according to the 
invention, again a triangular model of the pulse shape is 
used, and the actual pulse shape is a square-root -raised 
cosine pulse shape as depicted in figure 6. 

In this embodiment, however, the position estimate of the 
pulse peak is first calculated entirely as described with 
reference to figure l. A compensation of the difference 
of shapes 61, 65 of the model and the real pulse is only 
carried out in a subsequent step. 

For this compensation, the estimated position Xp is 
modified based on the amplitudes of the samples on both 
sides of the peak 62 of the pulse 65, and on the shape of 
the pulse 65. For instance, when the method described 
with reference to figure 1 is applied to the example of 
figure 6, it results in a negative error. This error can 
be reduced by a multiplication factor that is produced by 
the amplitudes of the samples at positions x 0 , Xi and x 2 . 

It is to be noted that the described embodiment 
constitutes only one of a variety of possible embodiments 
of the invention. 
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Claims 



1. Method for determining the position (xp) of a peak of 
a pulse in a signal received at a receiver, said 
method comprising: 

- taking samples of said received signal; 

- determining at least three samples, of which at 
least one has a signal strength exceeding a 
threshold value; and 

- determining the position (Xp) of said pulse peak 
based on an interpolation of at least two of said 
determined samples, which at least two samples are 
selected based on the signal strengths of said at 
least three determined samples, and which 
interpolation includes an evaluation of the signal 
strength of said at least two samples. 

2. Method according to claim 1, wherein different 
equations for said interpolation are provided for 
different distributions of the signal strengths of 
said at least three determined samples. 

3. Method according to claim 1 or 2, wherein said at 
least two samples are selected based in addition on a 
model for said pulse shape (31,41,51,61). 

4. Method according to one of the preceding claims, 
wherein equations for said interpolation are 
determined based on a model for said pulse shape 
(31,41,51,61) . 
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Method according to claim 3 or 4, wherein said model 
of said pulse shape (31,41,51,61) has a triangular 
shape . 



6. 



Method according to claim 5, wherein, in case the 
signal strength A(xo) of a determined sample at a 
first position x 0 is smaller than the signal strength 
A(xi) of a determined second sample at a following 
second position x x and the signal strength A(xi) of 
said second sample is larger than the signal strength 
A(x 2 ) of a determined third sample at a following 
third position x 2/ the position of said peak is 
estimated to be: 



Method according to claim S or 6, wherein said pulse 
has a width of about two chips, wherein said sampling 
rate is two samples per chip, and wherein, in case 
the signal strength A(x 0 ) of a determined sample at a 
first position Xo is smaller than the signal strength 
A(x x ) of a determined second sample at a following 
second position x x and the signal strength A(x x ) of 
said second sample is smaller than the signal 
strength A(x 2 ) of a determined third sample at a 
following third position x 2 , the position of said peak 
is estimated to be: 



4*,) 1 . . 

4O^ 0 )J chlpe - 



8. Method according to one of claims 6 to 7, wherein 
said pulse has a width of about two chips, wherein 
said sampling rate is two samples per chip, and 
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wherein, in case the signal strength A(x 0 ) of a 
determined sample at a first position x 0 is larger 
than the signal strength A(x x ) of a determined second 
sample at a following second position x x , the position 
of said peak is estimated to be: 



wherein Afx^) is the signal strength of a determined 
third sample preceding said first sample at a third 
position x-i, 

9- Method according to one of claims 4 to 8, wherein a 
weighting of the signal strengths of samples used in 
said interpolation is performed before said 
interpolation based on the known deviations between 
said model of said pulse shape (61) and the real 
pulse shape (65) . 

10. Method according to one of claims 4 to 8, wherein a 
correction of a position (x p ) determined based on said 
interpolation is performed based on the known 
deviations between said model of said pulse shape 

(61) and the real pulse shape (65) and based on the 
signal strengths of said samples. 

11. Method according to one of the preceding claims, 
wherein said at least three samples are consecutive 
samples . 

12. Device comprising means for determining the position 
(Xp) of a peak of a pulse in a signal received at a 
receiver according to one of the preceding claims > 




02/12 '02 LUN 16:40 [N° TX/RX 6019] 



2. DEZ. 2002 16:44 COHAU3Z & FLORACK 



NR. 050 6. 26 

PCT/IB02/05059 



- 19 - 

13. Device according to claim 12, wherein said device is 
given by said receiver. 

14. Device according to claim 12, wherein said device is 
a device external to said receiver and comprises 
further means for receiving from said receiver 
information on said received signal. 

15. Device according to claim 14, wherein said device is 
a network element of a cellular communication system. 

16. Cellular communication system comprising a device 
according to one of claims 12 to IS. 
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Abstract 



The invention relates to a method for determining the 
position Xp of a peak of a pulse in a signal received at a 
receiver- In order to improve the accuracy of this 
determination, the method comprises taking samples of 
said received signal/ determining at least three samples, 
at least one of which has a signal strength exceeding a 
threshold value, and determining the position Xp of the 
pulse peak based on an interpolation of at least two of 
the determined samples, which at least two samples are 
selected based on the signal strengths of the at least 
three determined samples, and which interpolation 
includes an evaluation of the signal strength of the at 
least two samples . The invention relates equally to a 
device and to a cellular communication system comprising 
means for realizing this method. 



For publication: Figure 1 
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take samples of received impulse 
response 



select sample at position xb which 
exceeds predetermined threshold 
and has a signal strength A(xo) 



select two subsequent samples at 
positions x 1t X2, which have signal 
strengths A(xi), A(x 2 ) 



yes 




peak position = 


1 

X, + — 


A(x 2 )~ 4*o ) 




1 2 


A{x 2 )+A(x 0 ) 





peak position = 



A( Xl )~A{xJ 
2lA{x x ) + A{x_ x )_ 



peak position = 


1 


AM 




x 0+I 


_A(x l )+A(x 0 ) m 





FIG. 1 
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FIG. 3 
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Estimation 




FIG. 4 




FIG. 5 



02/12 '02 LUN 16:40 [N° TI/RX 6019] 



, , 2. DEZ. 2002 16:44 



COHAUSZ & FLORACK 



4/4 



NR. 850 S. 30 

PCT/IB02/05059 




02/12 "02 LUN 16:40 [N° TX/RX 6019] 



